Monte Carlo Simulation for Elastic Energy Loss of 
High-Energy Partons in Quark-Gluon Plasma 



J. Auvinen , K. J. Eskola, H. Holopainen and T. Renk 

Department of Physics, P.O. Box 35, FI-40014 University of Jyvaskyla, Finland 
Helsinki Institute of Physics, P.O. Box 64, FI-00014 University of Helsinki, Finland 

E-mail: jussi . a.m. auvinenQjyu. f i 

Abstract. We examine the significance of 2 — > 2 partonic collisions as the suppression 
mechanism of high-energy partons in the strongly interacting medium formed in 
ultrarelativistic heavy ion collisions. For this purpose, we have developed a Monte 
Carlo simulation describing the interactions of perturbatively produced, non-eikonally 
propagating high-energy partons with the quarks and gluons from the expanding QCD 
medium. The partonic collision rates are computed in leading-order perturbative 
QCD (pQCD), while three different hydrodynamical scenarios are used to model the 
medium. We compare our results with the suppression observed in ^/snn = 200 GeV 
Au-f-Au collisions at the BNL-RHIC. We find the incoherent nature of elastic energy 
loss incompatible with the measured data and the effect of the initial state fluctuations 
small. 

1. Introduction: A substantial suppression of high-energy hadrons has been measured 
in = 200 GeV Au-hAu coUisions at the BNL-RHIC p. This is beheved to be a 
consequence of the energy loss of hard partons traversing a strongly interacting medium. 

To study the relevant physics of this phenomenon in as detailed manner as possible, 
we have developed a Monte Carlo (MC) simulation for the hard parton's interaction with 
the medium [2 [3], H]. This is similar to the (perhaps even more ambitious) MC mod- 
els JEWEL (Jet Evolution With Energy Loss) [S| and MARTINI (Modular Algorithm 
for Relativistic Treatment of heavy loN Interactions) [6]. However, while JEWEL and 
MARTINI include both elastic and radiative energy-loss components, we concentrate 
purely on the elastic energy loss for the time being. 

2. The model: We model the elastic energy loss of a hard parton by incoherent partonic 
2 — 7- 2 processes in pQCD, with scattering partners sampled from the medium. Three 
different hydrodynamical scenarios are used to model the QCD medium: i) a. (1-1-1)- 
dimensional hydro [7] with initial conditions from the EKRT model [8] for central heavy 
ion collisions, ii) a. (2-1-1 )-dimensional hydro [9] with a smooth sWN profile obtained 
from the optical Glauber model for non-central collisions, and Hi) an event-by-event 
hydro [9] with an eBC profile |10] from the Monte Carlo Glauber model to study the 
effects of the initial state density fluctuations. 
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Our approach is based on the scattering rate Vi{pi,u{x),T{x)) for a high-energy 
parton of a type i with 4- momentum pi, accounting for all possible partonic processes 
ij — )■ kl. The flow 4-velocity u{x) and the temperature T{x) of the medium are given 
by the hydrodynamical model. In the local rest-frame of the fluid, we can express the 
scattering rate for a process ij — )■ kl as follows [2j : 

\ roo riE\E2 

Tij^ki = ^Q^2Ef y_^^^2/i(^2,T) J^^^ ds[saij^ki{s,m% (1) 

Here Ei is the energy of the high-energy parton i in this frame and E2 is the energy 
of the thermal particle j with a distribution function fj{E2,T), which is the Bose- 
Einstein distribution for gluons and the Fermi-Dirac distribution for quarks. The 
scattering cross section aij^ki{s,m'^) depends on the standard Mandelstam variable s. 
A thermal-mass-like overall cut-off scale m = SmQsT is introduced in order to regularize 
the singularities appearing in the cross section when the momentum exchange between 
partons approaches zero. Here Qs is the strong coupling constant, which we keep fixed 
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with momentum scale. The free parameters of our model are thus Sm and 

The hard parton is propagated through the plasma in small time steps At. At each 
step, the probability for a collision is given by the Poisson distribution 1 — e"'"'^*. We 
always assume there is no significant interaction between the high-energy parton and the 
fully hadronic medium, and thus no collisions happen in regions with temperature below 
the decoupling temperature T^ec- The medium- modified distribution of high-energy par- 
tons obtained in the end can be convoluted with a fragmentation function to calculate 
the nuclear modification factor Raa{Pt, y, (f>) = ^N^J^^^liX/^idyd,^ • 

3. Results: In the following, our interest is in the high-PT neutral pions produced in 
^/sNN = 200 GeV Au+Au collisions. To achieve roughly the right amount of nuclear 
modification in the 0-10% centrality bin and to emulate also the incoherent higher-order 
processes, we set = 0.5 and Sm = 1-0. 

The simulation results for the 0-10%, 40-50% and 50-60% centrality bins, compared 
with the measurements, are shown in Fig. [H While the P-r-behaviour of the obtained 
Raa is compatible with the data within the studied transverse momentum range, it 
is clear from the figure that our model cannot reproduce the reaction plane angle 
dependence seen in the PHENIX experiment. Also the inclusive, angle- averaged 
nuclear modification factor fails to match with the experimental data: The computed 
suppression decreases too slowly as one advances to the more peripheral collisions. 

In Figure |2] we compare the angular dependence of the partonic Raa for fluctuating 
initial state geometry with the result for smooth initial conditions. While the variation 
in Raa between events can be considered notable, the angular variation within a sin- 
gle event is rather weak. The average over 20 events with fluctuating initial conditions, 
keeping the same value = 0.5 for the fluctuating and smooth cases, equals the smooth 
initial condition scenario with fairly good accuracy. In non-central collisions case the 
average of the events with fluctuating initial conditions is found to be systematically 
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Figure 1. (Color online) Left panel: The 7r° nuclear modification factor for 0-10% 
centrality, averaged over the reaction plane angle. Middle and right panel: The tt" 
nuclear modification factor dependence on the reaction plane angle A</) for 40-50% 
(middle panel) and 50-60% centrality (right panel). The simulation points (solid 
squares and triangles) are connected with lines to guide the eye. The PHENIX data 
are from [11] (0-10% centrality, open squares) and [12] (40-50% and 50-60% centrality, 
open triangles). Colored bars with small cap represent statistical errors; black bars 
with wide cap are systematic errors. 



above the smooth initial conditions curve, however. 

4- Summary: Our result for Raa{4>) in non-central collisions demonstrates that a purely 
incoherent energy-loss framework contradicts the present RHIC data. The weak sen- 
sitivity of the elastic energy loss model to the angle- dependent observables is clearly 
seen also in the fluctuating initial state study. In the central collisions, no difference is 
seen between the fluctuating and the smooth initial conditions when an average over 20 
events has been taken. In the non-central collisions the fluctuating conditions do appear 
to produce somewhat smaller suppression compared to the smooth background. 
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Figure 2. The partonic nuclear suppression factor Raa for central 200 AGeV Au-Au 
collisions at px — 10 GeV as a function of the angle of outgoing partons with respect to 
the event plane. Left panel includes plots for smooth initial conditions, for four different 
events with fluctuating initial conditions and for an average over 20 fluctuation events 
in 0-10% centrality. Nuclear suppression factor for tt*^ at Pt = 10 GeV in one event 
is also displayed. Right panel includes plots for smooth initial conditions and for an 
average over 20 fluctuation events in 30-40% centrality. 
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